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ABSTRACT 

We have used the Spitzer Space Telescope Infrared Spectrograph (IRS) to 
observe the 5-37 /im thermal emission of comet 73P/Schwassmann-Wachmann 3 
(SW3), components B and C. We obtained low spectral resolution (R~100) data 
over the entire wavelength interval, along with images at 16 and 22 yum. These 
observations provided an unprecedented opportunity to study nearly pristine 
material from the surface and what was until recently the interior of an ecliptic 
comet - cometary surface having experienced only two prior perihelion passages, 
and including material that was totally fresh. The spectra were modeled using a 
variety of mineral types including both amorphous and crystalline components. 
We find that the degree of silicate crystallinity, ~35%, is somewhat lower than 
most other comets with strong emission features, while its abundance of amor- 
phous carbon is higher. Both suggest that SW3 is among the most chemically 
primitive solar system objects yet studied in detail, and that it formed earlier or 
farther from the sun than the bulk of the comets studied so far. The similar dust 
compositions of the two fragments suggests that these are not mineralogically 
heterogeneous, but rather uniform throughout their volumes. The best-fit par- 
ticle size distribution for SW3B has a form dn/da~a^^'^, close to that expected 
for dust in collisional equilibrium, while that for SW3C has dn/da~a^'^°, as seen 
mostly in active comets with strong directed jets such as C/1995 01 Hale-Bopp. 
The total mass of dust in the comae plus nearby tail, extrapolated from to the 
field of view of the IRS peakup image arrays, is 3-5 x 10® kg for B and 7-9 x 10® 



^ Guest observer, Spitzer Space Telescope 

^ Guest Observer, Infrared Telescope Facility 

•^Also Department of Physics, University of Cincinnati 
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kg for C Atomic abundances derived from the spectral models indicates a deple- 
tion of O compared to solar photospheric values, despite the inclusion of water 
ice and gas in the models. Atomic C may be solar or slightly sub-solar, but its 
abundance is complicated by the potential contribution of spectrally featureless 
mineral species to the portion of the spectra most sensitive to the derication of 
the C abundance. We find a relatively high bolometric albedo, ~0.13 for the 
dust, considering the large amount of dark carbonaceous material, but consistent 
with the presence of abundant small particles and strong emission features. 

Subject headings: comets: general comets: individual (73P/Schwassmann-Wachmann 
3) methods: data analysis - techniques: spectroscopic 



1. Introduction 

A major goal of planetary science is to develop a theoretical framework that is capable 
of reproducing the characteristics of our solar system today. Comets form a direct link to 
the physical and chemical processing that occurred when our solar system was was only a 
few Myr old (0.02% of its present age). These primitive materials consist of ices of abundant 
volatile molecules (water, methane, ammonia, carbon monoxide, carbon dioxide), and dust - 
highly refractory inorganic materials (rock-forming minerals), as well as organic materials of 
generally intermediate volatility. Of particular interest are the mineralogy, size distribution, 
and degree of crystallinity of the dust. These provide key information about the condensation 
of mineral types from the young solar nebula, the dust-gas chemistry, as well as possible 
large-scale transport within the nebula. 



1.1. Origin of Comets and their Dust Material 

Comets fall into two broad classes, based on their orbital dynamics. The nearly isotropic 
comets (NICs, previously called dynamically-new comets and Oort cloud comets) are those 
entering the inner solar system from nearly isotropic directions on highly eccentric orbits. 
They are thought to have been formed and ejected from the Jupiter-Uranus region, and 
are currently stored in the Oort cloud. The ecliptic comets (ECs), by contrast, are those 
with low orbital eccentricities and with inclinations close to the ecliptic, and which include 
objects significantly affected by the gravity of Jupiter, so that they are in resonant or near 
resonant orbits with the giant planet, the Jupiter-family comets (JFCs). They are commonly 
believed to have formed further away from the Sun than the NICs were, near the outer edge 
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of the nebula/disk, and they might be expected to possess distinctly different compositions 
than the Oort cloud objects. Dynamical studies suggest that the ECs as well as many of 
the intermediate Halley-family comets (periodic non-ECs) and possibly some NICs come 
from the trans-Neptu nian scattered disk, a population still exhibiting these eccentricity and 



inclination patterns ([Duncan fc LevisonI 119971 : iFernandez et al.l l2004l : iLevison et al.l 12006 



Volk fc Malhotrall2008l ). Simulations of the orbital evolution of the pla nets indicate that sig - 
nificant migration (and possibly even positional swapping by planets (JTsiganis et al.l 120051 )) 
may have occurred during the time of comet formation, so that the distinction between the 
location of formation of the two comet populations may be blurred. 

The composition of the minerals that make up cometary dust will be a mixture of 
pre-solar grains that are unchanged remnants of interstellar medium (ISM), material that 
accreted into the molecular cloud that the Sun formed in, materials that condensed within 
the cloud and proto-solar nebula itself, and material which has undergone a sequence of 
condensation and evaporation and annealing in the circulating proto-planetary disk formed 
after the nebula collapsed to form the young Sun and its accretion disk. The material that 
condenses (or re-condenses) in the nebula will be subject to a complex set of gas-grain (and 
gas-gas) int eractions w hich can produce a radial stratification in the condensation sequence 
of minerals (JGaillll998l ). Because of this, it is possible that comets forming at different radial 
distances may have inherited different mineralogies. 

It is known tha t the grains in the interste llar medium have a very low degree of crys- 
tallinity, only ~1% (JKemper et al.l l2004j . |2005| ) . This contrasts to the much h i gher values 
observed in comets - exceeding 50% in some cases (JHarker et al.l |2002| . |2004J : iLisse et al. 



20071 ). Where, when, and how this transformation occurs is still debated, but annealing and 
sublimation/condensation in the hottest part of the solar nebula is implicated. Spatially 
resolved spectroscopy of the silicate band of dust in young stellar disk systems indicates that 
the inner few AU have larger and more crystalline silicates than the more distant regions 
(Ivan Boekel et al.l |2003[ ) . suggesting production in the inner regions and radial transport 
outward. 

A related aspect of the formation and processing of dust in the solar nebula is the 
availability of elemental species. Even the most "primitive" meteorites, the CI carbona- 
ceou s chondrites, are dep leted in carbon and oxygen compared to solar photospheric val- 
ues (JLodders et al.ll2009al Jbl). A similar pattern is seen i n the mass spectrorneter data on 
IP/Halley obtained with the Gio t to and VEGA space craft (JLangevin et al.lll987l : IJessberger et al 
19881 : IJessberger fc Kisselll Il99ll : iLawler et al.l Il989l ) , which presumably sampled material 
condensed further from the Sun than the meteorites. These are critical elements for the 
formation of silicates, ices, and organics materials in the solar system. A radial gradient 
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in the depletion pattern might provide clues about the condensation processes in the early 
solar nebula, and signal a lack of significant radial mixing. 

If there is a radial gradient in the characteristics of grain mineralogy and degree of 
crystallinity of the material in the early solar nebula, it would be manifested in differences 
between the grain characteristics of these two comet populations. Furthermore, if the grain 
characteristics of the comet-forming regions evolves while the cometesimals are growing and 
assembling, these bodies may be internally heterogeneous and possibly stratified. These 
differences should manifest themselves in differences among the grain characteristics of the 
comet populations. However, after many passages through the inner solar system, the sur- 
faces of EC comets will become de-volatilized, processed, and may no longer represent the 
material they formed from. Furthermore, if the mineralogy and degree of crystallinity in the 
comet-building region was evolving as the comet nuclei were being assembled, then detection 
of inhomogeneities within the nucleus would provide a window on the temporal changes in 
the comet-building zone. Any stratification in grain properties would be difficult to detect 
from a mere snapshot study of its current (and processed) surface. 

For these reasons it is important to sample the inside of comet nuclei, particularly those 
of the ECs. One method is to artificially excavate beneath the surface layers, as was done 
by the Deep Impact (DI) mission to comet 9P/Tempel 1 (hereafter Tl). Another approach 
is to study comets that have recently fragmented, exposing their more pristine interiors for 
the first time. 

During the past decade or so, a number of ECs have been observed to fragment: 
D/1993 Shoemaker-Levy 9, 16P/Brooks 2, 51P/Harrington, 57P/du Toit-Neujmin-Delporte, 
73P/Schwassmann-Wachmann 3 (hereafter SW3), lOlP/Chernyk h, 102P /Shoemak er-Holt 1, 



174P/Echeclus, 205P/Giacobini, and P/2004 V5 (LINEAR-Hill) flFernandezll2009f ). Because 
disrup ted comets may be res ponsible for the production of most of the dust in the zodiacal 
cloud ( Nesvorny et al.ll2010l ). understanding the physical and chemical properties of these 



comets is important. As of 2006 there were no mid-IR spectra of a fragmented EC suitable 
for the necessary analysis. However SW3 survived a fragmentation event in 1995, and two 
orbits later was poised for a close approach to the Earth. This provided an unprecedented 
opportunity to study one of these objects in detail. 

In this paper, we describe observations of the B and C fragments of SW3 using the In- 
frared Spectrograph (IRS) of the Spitzer Space Telescope. Although Spitzer was more distant 
from the comet than the Earth at its closest approach, its superior sensitivity, wavelength 
coverage, and lack of telluric spectral contamination provided us with the best opportunity 
to examine this important object in the mid-IR. Additional supporting ground-based visual, 
near-IR, and mid-IR observations that were useful for interpreting the Spitzer data are de- 
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scribed in the Appendix. We will also compare the results of the analysis of the Spi tzer 



observations with other independent mid-IR observations of SW3 (iHarker et al.l |2010| ) as 
well as other comets studied spectroscopically and using in situ sampling techniques. 

In this study we will be able to: 

• obtain information on the mineral content and degree of crystallinity on (potentially) 
the most pristine cometary material yet measured through remote sensing 

• compare the abundances of the various minerals between two cometary fragments that 
represent different fractions of material formerly buried below the surface of the nucleus 

• compare the atomic abundances of SW3 with other comets studies both spectroscop- 
ically and in sutu, primitive meteorites, and the Sun, to investigate the pattern of elemental 
depletions among these solar system objects 



2. Observations & Spectral Processing 

IRS observations of SW3-C and SW3-B were obtained on 17 March 2006 UT, and 17 
April 2006 UT, respectively. For both objects, image maps were obtained using the Blue 
Peak-Up array (A = 13.3 — 18.7/im, hereafter referred to as the 16 /xm array) and Red 
Peak-Up array (A = 18.5 — 26.0/im, hereafter the 22 /zm array). The primary spectroscopic 
observations were obtained in the "stare" mode with the low spectral resolution gratings. 
The various slits, spectral orders, wavelength coverage, and spectral resolutions R = A/AA 
of the IRS are listed in Table 1. In each case, the object was acquired with the 22 /xm array, 
then offset to the IRS slit for the spectral observations. 

SW3-C was observed at a hehocentric distance r/i=1.47 AU, and distance A= 0.78 AU 
from Spitzer. The 16 & 22 /im image maps indicated that this component consisted of a 
single fragment, with no evidence of smaller fragments within the field of view. 

SW3-B was observed at r/i=1.19 AU, and a distance A= 0.44 AU from Spitzer. The 
observations were o btained within a day of the first set of Hubble Space Tele scope images 
( Weaver et al.ll2006l) . when the tail of SW3-B consisted of many sub-fragmentq^. Ilshiguro et al 



( I2OO9I ) reported detecting over 150 fragments using the Subaru Suprime camera, in addition 



^The 16 /^iTL Spitzer IRS images have a point spread function (PSF) of 3.6 arcsec, which at the distance 
of SW3-B (0.44 AU) was too broad to detect the individual "crumbs" seen in the Hubble images, which had 
superior spatial resolution due to its larger collecting mirror, shorter wavelength used (0.6 /ini), and closer 
distance at the time of observation (0.22 AU). 
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to the main int act body. The frag mentation that produced these pieces probably occurred 



around 1 April ( IGreen et al.ll2006[ ). about 16 days before the Spitzer observations. 



Observations of SW3-G were also attempted, but failed to produce any data of sufficient 
quality for extraction and analysis. Hubble images obtained less than 30 hours after the failed 
Spitzer acquisition revealed G to have consisted entirely of a swarm of sub-fragments that 
were apparen tly the result of a m aior disintegration event a few days prior to the Hubble 



observations (IWeaver et al.l 120061 ). The failure was thus likely due to G being too diffuse, 



too faint, or both. 

In addition to the low-resolution "stare" observations, observations were also made in 
the low-resolution "spectral mapping" mode to insure detection (albeit at a reduced signal- 
to-noise ratio) if the ephemeris was in error g, to obtain spatial-spectral information on the 
main components, and possibly measure the spectra of newly-shed fragments. A subset of 
the spectral maps were used to determine the slit loss correction factors for extracting the 
spectra from the "stare" observations (see Sec. 2.2). High spectral resolution measurements 
were also made in order to recover spectral information should the low-resolution spectra be 
saturated. In this paper we will concentrate solely on the low-resolution "stare" observation^^' 



2.1. Image Maps 

As part of the standard target acquisition process for Spitzer's IRS, images are obtained 
from the peak- up arrays, are then used to center the target on the spectrograph slit. For 
SW3, these were saturated in the comae and not used for extracting science information. 
For this work, we instead used image maps constructed from a number of additional, short 
6 sec exposure, dithered images we requested in each of the arrays. 

The IRS peak-up array images were bad pixel masked, re-projected onto the tan- 
gent plane, and rnosaic ked in the rest frame of the comet using the MOPEX software 



( iMakavoz &: Kahnll2005[ ) version 18.2.0. The 22 yum image maps covered clean backgrounds, 
but the 16 /im images map did not. We estimate the 16 /im background by measuring the 
faintest edge of the images and multiplying by a factor of 0.97. These edges are in the direc- 



^ Observations needed to be planned weeks before the initial recovery of the various fragments of the 
comet. 

•^The high resolution spectra do not include data shortward of 9.9 /im, are noisier than the low resolution 
data, and (as of the time of the analysis) suffered severe order-to-order mis-matches for extended sources 
that made their merger problematic. However, the extended source calibrations for the high-resolution data 
are improving and these data might be tractable in the near future. 



tion of the Sun, i.e., they do not coincide with the comet tails. The factor of 0.97 was derived 
by comparing the faint, near- nucleus image edge of the 22 /im images to their background 
observations. This method is likely accurate to 0.1 MJy sr~^ pixel"^. The 16 fim mosaics of 
fragments B and C are presented in Fig. [H 



2.2. Extended Source Calibration and Spectral Extraction 

IRS spectra are calibrated with spectroscopic observations of point sources. The IRS 
slits are narrow with respect to the PSF of the telescope, i.e., they do not encompass 100% of 
the PSF at any wavelength and the fraction of the PSF encompassed varies with wavelength. 
Thus, the low-resolution IRS slits block up to 35% of the flux of a point source at the longest 
wavelengths, but as little as 17% at the shortest wavelengths. Because the IRS instrument is 
calibrated with observations of unresolved stars, the slit-losses may be ignored when working 
with calibrated observations of point sources. Ignoring the slit-loss does, however, affect the 
spectral shape of calibrated observations of sources larger than the PSF of the telescope. In 
the case of an infinitely extended source the slit-losses are zero, but the calibration assumes 
losses equivalent to a point source. Therefore the resulting "calibrated" spectrum will not 
accurately represent the spectral shape of the infinitely extended source. 

The thermal emission from comets is comprised of two parts: 1) emission from the 
nucleus, generally spatially unresolved, and 2) emission from the dust, which is a region of 
radially and azimuthally varying surface brightness. To properly flux-calibrate our comet 
spectra, we compute the slit-losses for the observed surface brightness distribution, and the 
orientations and sizes of each IRS slit. No correction was applied for the presence of the 
nuclei of either fragment, as the emission from the nuclear surfaces was negligible compared 
to the emission from the huge surface areas of fine coma dust, (see Sec. 3.1). 

We began by measuring the spatial profile perpendicular to the observed tail of each 
comet in the 5-7 /x m region of the IR S data cubes built from the spectral maps using the 



CUBISM software (ISmith et al.l 120071 ). The 6 fxm. region of the data cubes was chosen for 
its high spatial resolution and good quality signal. We fit photometric cuts perpendicular to 
the tail with both Moffat and Gaussian functions, but the Moffat functions fit the images 
better. From the best-fit Moffat parameters we generated a model 6 fim image of the comet 
at an image scale of 0.1835 arcsec pixel"^, i.e., one-tenth the image scale of the IRS peak- 
up arrays. The super-resolution model image is convolved with the synthetic IRS 16 /im 
peak-up array PSF from a 300 K blackbody spectrum, to roughly approximate the the 
spectral shape of the thermal emission. This PSF was generated at an image scale of 0.1835 
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arcsec pixel"^ with version 2.0 of the Tiny Tmv/ Spitzer softwarqj. With the addition of a 
point source component for fragment C, the surface brightness distribution of the convolved 
super-resolution model images approximately match the peak-up imaging observations. The 
parameters of the Moffat fits, plus the additional point source for fragment C, form the basis 
of our model comets with which we derive the slit-losses of each comet. 

We measured the slit-losses of the low resolution modules and the short wavelength, high 
resolution module by aligning a rectangular aperture on each comet's nucleus, reproducing 
the orientation and size of each slit with respect to the comet. The flux contained within the 
slit is computed for A) the model image, and B) the model image convolved with a synthetic 
300 K blackbody PSF computed for the IRS instrument at wavelength A, and smoothed 
with a 1.1 (SL) or 1.3 (LL) pixel wide boxcar function. Without the smoothing, a Tiny Tim 
3000 K PSF would not match the observed width of calibration stars processed by CUBISM. 
The slit-loss is then 1 — B{\i)IA. Figs. |2] & [3] presents our model images, the orientations 
of the IRS slits, our computed slit-losses, and the slit-losses for a point source as provided 
by the Spitzer Science Center in the Spitzer IRS Custom Extraction (SPICEJfl software's 
calibration flies. To produce the flnal flux calibrated spectra, we take the extended source 
calibrated spectra from SPICE (which assumes an inflnitely extended surface brightness) 
and divide it by B{Xi)/A. We followed this prescription for all but the SL observations of 
Fragment SW3-B, because the tail of fragment B is smoothly varying over many slitwidths, 
and the SL slit is nearly aligned perpendicular to the tail, we elected to use the SPICE 
extended source calibration for the SL spectrum of this fragment. 

The resulting spectra of SW3-B and SW3-C are shown in Figs. |l]&[5l For both objects 
the flux levels for the LL2 (14.0-21.3 /im) and SL2 (5.2-8.7 /xm) modes of the IRS have been 
normalized to that of SLl (7.4-14.5 /im) at the wavelengths where they overlap, while LLl 
(19.5-38.0 fim) has been normalized to the re-scaled flux of LL2. The overall difference in 
spectral shape between SW3-B and SW3-C occurs because they were observed at different 
heliocentric distances (1.19 AU for B and 1.47 AU for C) with different local insolation flux 
densities, and because the particle size distributions for the dust flowing out into the coma 
were different as well for the two comets - the B-fragment coma contained many more large 
particles, pres umably due to the p rocesses creating the large fragments streaming down the 
tail by Hubble ( Weaver et al.ll2008l ). Both result in producing somewhat different equilibrium 



temperatures for the grains. 



*Tiny Tim./ Spitzer is available from the Spitzer Science Center: 



http://ssc.spitzer.caltecli.edu/ 



^ http://ssc.spitzer.caltech.edu/ciataanalysistools/tools/spice/spiceusersguide/ 
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3. The Dust Content of Schwassmann-Wachmann 3 
3.1. Spectral Modeling 

The infrared emission from a collection of dust near a star is given by 

Fx,mod = ^^ Bx{Ti{a,rh))Qabs,M^^)'^(^'^ — j-^da (1) 

where T is the temperature for a particle of radius a and composition i at heliocentric 
distance r/j, at a distance A from the observer, Bx is the blackbody radiance at wavelength 
A, Qabs is the emissivity (emission efficiency) of the particle of composition i at wavelength 
A, dn/da is the differential particle size distribution (PSD) of the emitted dust, and the sum 
is over all species of material and all sizes of particles for the dust. Spectral analysis consists 
of calculating the emission flux for a model collection of dust, and comparing the calculated 
flux to the observed flux. The emitted flux depends on the composition (location of spectral 
features), particle size (feature to continuum contrast), and the particle temperature (relative 
strength of short vs. long wavelength features). 

3.1.1. Scattered Light and Nucleus Contamination 

Prior to the modeling, the effects due to reflected sunlight and nucleus thermal emission 
need to be checked, and if necessary, remove d. Using pho t ometr y obtained on 28 March with 



the Spitzer Infrared Array Camera (IRAC) iReach et al.l ( 120091 ) found the scattered light to 



be negligible at the 6% level at 4.5 /xm and dropping rapidly at longer (i.e. IRS) wavelengths; 
the same study also provided an accurate check of the Fi,(23.7 yum)/F,^ (7.7 microns) in our 
extracted IRS spectra. 

We independently checked the degree to which scattered light might affect the shortest 
wavelengths of the Spitzer IRS (5 fun) of component C with ground-based spectrophotome- 
try covering 0.44-13.5 /zm obtained 19 May using The Aerospace Corporation's Broad-band 
Array Spectrograph System (BASS) and its CCD guide camera on NASA's Infrared Tele- 
scope Facility (IRTF), and found it to be less than 11% and likely smaller (see Appendix). 
Because the main contributor to the 5 /im emission will be amorphous carbon (see below) its 
abundance might be too high by 6%-10%. Other species will be affected to a lesser extent. 



Barker et al.l ( 120101 ). using 10 /xm imaging obtained with the Michelle instrument on 



the Frederick C. Gillett Gemini North (hereafter Gemini North) telescope determined the 
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contribution of the nuclei in a 0.6" x 1.0" synthetic beam, and found it to be 1.4%±0.2% 
for B and no detectable contribution to C (both at a 99% confidence level). Our Spitzer 
spectra were obtained with a much wider slit and at over twice the distance to the comet, 
implying an even larger beam/slit dilution factor for an unresolved point source factor, and 
so we applied no correction for the nuclear emission in this work. 



3.1.2. Composition 

The spectral model used was the same one applied to the DI spectra of Tl and the 
Infrared Space Telescope ( ISO) Short Wavele ngth Spectrograph (SWS) observations of Hale- 
Bopp (hereafter HE) by iLisse et al.l (120071 ). To determine the mineral composition the 
observed infrared emission is compared with a linear sum of laboratory thermal infrared 
emission spectra. As-measured emission spectra of randomly oriented 1 /xm-sized powders 
(produced by grinding of representative mineral samples) are utilized to determine Qabs 
for a statistical collection of lowest free-energy dust particles. The list of species tested, 
motivated by their reported presence in interplanetary dust particles, meteo rites, in situ 
comet measurements, young stellar objects, and debris disks fJLisse et al.l 120071 ) includes fer- 
romagnesian silicates of various compositions (forsterite, fayalite, clino- & ortho-enstatite, 
augite, anorthite, bronzite, diopside, & ferrosilite); silicas (quartz, cristobalite, trydimite); 
phyllosilicates (such as saponite, serpentine, smectite, montmorillonite, & chlorite); sulfates 
(such as gypsum, ferrosulfate, & magnesium sulfate); oxides (including various aluminas, 
spinels, hibonite, magnetite, & hematite); Mg/Fe sulfides (including pyrrohtite, troilite, 
pyrite, & niningerite) ; carbonate minerals (including calcite, aragonite, dolomite, magne- 
site, & siderite); water ice (clean and with carbon dioxide, carbon monoxide, methane, & 
ammonia clath rates); carbon dioxide ice; graphitic and amorphous carbon; and neutral and 
ionized PAHs JPrain fc Lill2007l : llisse et al.lbood . l2007h . 



3.1.3. Particle Size and Dust Mass 

Particles of 0.1-2000 /im in radius are used in fitting the 5—37 /im Spitzer data, with 
particle size effects on the emissivity assumed to vary as 



1 - Qa,X = [1 - Q 



l^m,,A 



(a/l/xm) 



(2) 



The particle size distribution (PSD) is fit at approximately logarithmic steps in radius, 
i.e., at [0.1, 0.2, 0.5, 1, 2, 5,100, 200, 500, 1000, 2000] /xm. Particles of the smallest sizes 
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have emission spectra with very sharp features, and httle continuum emission; by contrast, 
particles of the largest sizes are optically thick, and emit only continuum radiation. Spectra 
with exceptionally strong spectral features thus have a significant abundance of small (~ 1 
/im) grains. 

As discussed in §2 above. Fragment 73P/SW3-C was observed at A= 0.78 AU, Vh =1.47 
AU on 17 March 2006. The best-fit PSD found from fitting the IRS spectrum, after division 
by a 280 K greybody to convert it into an emissivity spectrum (Figs. El&ED, is a relatively 
steep one similar to the one found for the very small particle dominated Hale-Bopp coma 
dust: dn/da = a~^'° (dn/dlog (m) = m~^°, 0.1 /im < a < 1 mm), but with an additional large 
particle component, a Gaussian (or wide delta function) a t ~ 1 mm to fit the observe d long 



wavelength continuum. A similar bump was derived by IVaubaillon fc ReachI (|2010| ) from 
Spitzer images of the comets and their associated meteoriod streano. The total in-beam 
surface area and mass is 38 ± 2 km^ and 1.6 ± 0.1 x 10^ kg, respectively (equivalent to an 
11.5 m radius body of 2.5 g cm"^ density). 

Fragment 73P/SW3-B was observed at A = 0.44 AU, Vh = 1.2 AU on 17-19 Apr 20ofl 
The best-fit PSD found from fitting the IRS spectrum, after division by a 300 K greybody to 
convert it into an emissivity spectrum (Figs. [8] &[9]), is a moderate one similar to those found 
for collisional equilibrium systems: dn/da = a~^'^ (dn/dlog (m) = m~°'^^, 0.1 yUm < a < 1 
mm), but with an additional large particle component, a Gaussian (or wide delta function) 
at ~1 mm to fit the observed long wavelength continuum. The total in-beam surface area 
and mass, dominated by contributions from the largest particles, is 12 ±1 km^ and 4.1 it 0.2 
X 10^ kg, respectively (equivalent to an 7.4 m radius body of 2.5 g cm~^ density), somewhat 
less than found in the Fragment C coma, despite the B-Fragment being 20% closer to the 
Sun at the time of observation. The PSDs of both B and C components are shown in Fig. [TOl 

Extrapolating the in-beam values to the rest of the comae using the peak up images 
yields a total coma mass of 4-7 x 10* kg and 7-9 x 10® kg for B and C, respectively, contained 
in the 80 x 56 arcsec field of view of the blue peakup image array. For both fragments, the 
equivalent si ngle body est i mates are probably lower limits. Radar measurements of both 



fragments by iHowell et al.l (120071 ) indicated the presence of cm-sized "gravel" in the comae, 
which would raise the mass of dust in the comae. For a dn/da ~ a"^*^ PSD, we estimate 
the mass will increase as log(amax), while for dn/da ~ a~^'^ PSD, the total coma mass will 
increase as (amax)^'^, where amax is the radius of the maximum sized particle. 



A simila r 1 mm bump was al so present in the PSD of the dust of IP/Halley ( McDonnelll Il99l[ ) and 



81P/Wild 2 (|Tuzzolino et al.ll2004l ). 



^including all spectral modes 
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The slope of the particle size distributions dn/da = a ^'^ and dn/da = a ^'^ for C and 



B respectively, are sin i ilar to the va.lue of dn/da = a 



1997 



-3.6 to 4.0 



found for HB (j Williams et al. 



Lisse et al.lll999l : iHarker et al.ll2002[ ). which had the most abundant population of small 



grains measured in any comet. The high abundance of small particles may explain the unusu- 
ally h igh polarization of SW3 compared to that of other ECs (JHadamcik fc Levasseur-Regourd 
20091 ). A small grain population usually correlates with a strong silicate emission band at 
10 /xm, and while that band strength in both B and C are stronger than that typically seen 
in other ECs, they are weaker than that of HB. For fragment C the PSD is at least as steep 
as that of HB, but a combination of a higher abundance of carbon and a greater fraction of 
its silicates in the amorphous phase (discussed below), will result in a more muted silicate 
band. For fragment B, with its shallower PSD than HB, the inclusion of a large number of 
"boulders" in the spectrograph slit is likely add a significant amount of grain emission with 
a weak silicate band. 

The difference in the power law exponent in the size distributions between the two frag- 
ments does not necessarily point to an intrinsic difference in the origin of the grains in the two 
fragments measured. As discussed above, the profuse shedding of "boulders/gravel/crumbs" 
by B insured a larger collection of bigger particles than for C Furthermore, th ere is consider- 
able e viden ce for ongoiiig grain fragm entation in the outflowing dust in comets. iTuzzolino et al 
( 120041 ) and lGreen et al.l (120041 12007| ) reported that the Stardust Flux Monitor measurements 
of comet Wild 2 contained fine-scale structure that was best understood by particle fragmen- 
tation. Simil ar structure events were detected in the dust outflows of Halley by the VEGA 
instruments (jOberd 120041 ). Observations of the material ejected from Tl immediately af- 
ter the DI event are also consisten t with extended de-aggregagtion of the ejected material 



(iLisse et al.ll2006l: Keller et a 



20071) . The polarimetric imaging data on both large fragments 



of SW3 itself by I Jones et al.l ( l2008l ) indicated grain f ragmentation occurr ed as the dust flowed 



away from the nucleus. The radar measurements of iHowell et al.l ( 120071 ) of SW3 may also be 
consistent with fragmentation occurring in the coma. The dust fragmentation will be very 
stochastic in nature, and the different character of the large-scale fragmentation of the nuclei 
of B compared to C will further complicate the eventual particle size distributions derived 
from the spectroscopy. 



3.1.4- Particle Temperature 



Dust particle temperature is determined at the same log steps in radius used to deter- 
mine the PSD. Particle temperature is function of particle composition and size at a given 
heliocentric distance. The highest temperature for the smallest particle of each species is free 
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to vary in our model, and is determined by the best-fit to the data; the largest, optically thick 
particles (2000 /xm) are set to the local thermodynamic equilibrium temperature (LTE), and 
the temperature of particles of intermediate sizes is interpolated at each PSD point between 
these extremes. Good agreement was found using this method for the temperature of the 
Tl dust ejecta andthat expected from LTE, where Tmaa(best fit) = 1.4 Tlte, and by com- 
pa rison of the Tm.n .m and T/^.tr temperatures found for the HB coma dust vs those estimated 
by lCrovisier et al.l (Il997| ). 



3.2. Overall Mineralogy 



The overall compositional mineralogy of the dust released from the two largest and 
brightest fragments of the 2006 apparition, 73P/SW3-C and 73P/SW3-B, is similar, as 
illustrated in FigJTT] Both show that is dust dominated by, in order of abundance, amorphous 
carbon, amorphous silicates of near-pyroxene composition, Mg-rich crystalline forsterite, Fe- 
metal sulfides, water ice and gas. The major difference between the best-fit models of the 
two fragments resides in the form of water: in fragment B the effective weighted surface area 
of water vapor is larger than that of water ice, while this ratio is reversed in fragment C 
The relative strengths of the two fragments' water emission lines above the dust continuum 
at 6 yum is apparent in Fig. [T21 PAHs, carbonates and crystalline pyroxenes, found in 
comets Tl, HB, and 17P/Holmes ( JLisse et al.l 120071 : iReach et al.ll2010l ). are only marginally 
detected, if at all. The strong similarity between the released dust from the two different 
fragments is striking, in that they have very different emitted particle size distributions and 
temporal histories, with the B-fragment going into long extended outburst periods of high 
activity during April-May 2006, and with the likelihood that the B and C fragments have 
different vertical structure contributions from the original pre-split nucleus. Nevertheless, 
the derived Spitzer dust compositions are identical within the errors of the modeling, arguing 
for homogeneity of the original nucleus, consistent with the narrow band optical photometry 



and near-IR spectrosco py of gases released from the two bodies (jSchleicher fc Bairl 12008 
Dello Russo et al.ll2007l ). 



3.3. Silicates 



The SW3 silicates are dominated by their amorphous component of near-pyroxene com- 
position, much more so than the dust ej ected by the DI ex periment from Tl or the abundant 
dust flowing through the coma of HB (JLisse et al.l 120071 ). Of the crystalline silicates, only 
Mg-rich olivine (forsterite) is found in abundance for the best-flt mineralogical models of the 
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dust released from SW3 B & C. 

The presence of crystalline silicate s above tha t seen in the interstellar medium (ISM) 
is a common feature of cometary dust. iLisse et al.l (|2006[ ). analyzing the dust ejected from 
Tl during the Deep Impact event, found that both olivines and pyroxe nes were more than 



70% crystalline. In HB, the fraction was lower, between 30% and 60% (JHarker et al.l 12002 



Lisse et al.l 120071 ) In SW3, despite being much more primitive than HB or Tl (Sec. 4. 2), 
the dust is still far more crystalline than ISM silicates. The fractional crystallinity of the 
silicates in both components of SW3 derived from the Spitzer dat a is 25%-30% . An in depen- 
dent estimate of the crystallinity of the sihcate dust of SW3 by iHarker et al.l (120 lOf ). using 
spectroscopic observations of components B and C with the Mic helle mid-infrared sp ectro- 
graph on the Gemini North telescope, and the modeling code of iHarker et al.l ( 12002| ) gives 
the values of 33.5+^50% for B and 25.7+r^%. 



'-25.3 



4.3' 



The 8-13 yum silicate feature strength, 50% above the nearby continuum, is significantly 
higher than the 10-20% typically seen for ECs, but lower than that seen in HB or the post- 
impact spectrum of Tl. The change in band strength in the pre- and post-impact spectra of 
Tl is most likely caused by de-aggregation of large, optically thick (> 10 /im) porou s dust 
particles into their sub-micron to micron constituent components ( ILisse et al.l 120071 ). The 
same process might account for the strong (compared to other ECs) silicate feature in the 
SW3 fragments, and is presumably due to the same mechanism causing the fragmentation of 
SW3. The even greater strength of the silicate feature in HB and the post-impact spectrum 
of Tl is likely be related to the more energetic ejection mechanisms involved: the high-energy 
of the DI impactor and the violent outflow from HB compared to the more gentle disruption 
process in SW3. 



3.4. Amorphous Carbon 

Another striking result of the SW3 compositional modeling is the very large amount of 
amorphous carbon in the mix. The relative surface area for this component is huge, as can be 
seen in the very flat emissivity spectrum in the 5-9 /xm range. An i ndependent consiste ncy 
check of this result was found by applying the Mie dust models of IHarker et al.l (120041 ) to 
the Spitzer data, with similar results. There is much more C in the solid phase than for Tl 
or HB. 



For comparison, the most carbon rich meteorite, Tagish Lake, is 3.6±0.2 % (JBrown et al. 
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2000) carbon by masqj. Our results for the outflowing SW3 B & C dust is about an order 
of magnitude more C-rich by mass. It must be noted that as the emissivity spectra for 
amorphous C and nano-phase Fe are both very similar and featureless in the mid-IR, it is 
near-impossible to spectrally distinguish between the two and thus the total carbon abun- 
dance could be somewhat lower than our estimates. We can, however, put a weak upper 
limit on the nano-phase Fe content, assuming the elemental abundance of Fe to be less than 
or equal to solar (Figs. [13] and [H]) . This implies that nano-phase Fe is able to account for at 
most ~l/3 of the superabundant amorphous C signal. Additional emission from other spec- 
trally featureless materials, such as very large grains, might also be driving the C abundance 
high in the Spitzer spectra. We note, however, that the Giotto and VEGA mass spectrom- 
eter measurements of Halley have a lower carbon content (compared to the silicate-forming 
elements Si and Mg) than SW3, and are comparable or slightly higher than Tl, HB, or 
Holmes. 



3.5. Water Ice and Gas 

Even without application of compositional modeling, the stronger water gas lines in the 
Fragment B spectrum are clear in Fig. [TTJ This finding is somewhat s urprising, given the 



HST imaging of multiple icy bodies moving down the tail of Fragment B fjWeaver et al.ll2006 



20081 ). which would seem to imply abundant amounts of solid water ice. On the other hand, 
the distance from the Sun, 1.2-1.4 AU, and the gross overall best fit dust temperatures, 
280-300 K, both suggest that the majority of emitted water should be in the gas phase, and 
it is possible we are seeing the result of an outburst of water ice some hours or days after it 
has had a chance to warm up and evaporate. 

Taking note of the much finer size distribution (dn/da ~ a~^°) of emitte d dust for 



Fragm ent C, its identification as the nominal main nucleus (JLamy et al.l |2004J : ISekanina 



20051 ). providing all the other fragments, and its near-normal emission behavior (in terms 
of the temporal development of its light curve, and the morphology of its coma), the most 
plausible explanation we can find is that the water emitted from a well mantled (due to 
fallback of heavy dust fragments) Fragment C is supplied predominantly from subsurface fine 
ice grains which are slowly subliming as they fiy, while the water coming out of Fragment 
B is fiowing out mainly as water vapor, either from the fresh surface of the main fragment 
B body or the large optically thick (i.e. spectrally undetected) sub-fragments moving down 



^One sample had an abundance of 5.81 %, higher than any other known chondrite (JGradv et ahl 120021) . 
but this is still less than what we derived for SW3. 
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its tail. Evidence for abundant icy dust in the fragment B tail is provided b y the narrow 
pinch ing" morphological effect seen in the images of B seen in Fig. [T](see also 



larker et al. 



( l2010l )). Such morphology is predicted by the dynamical models of iLienI (jl996[ ). where the 
angular width of the observed tail is substantially and atypically narrower than that of the 
coma, due to ongoing ice sublimation. 

Finally, these interpretations are complicated by the fact that th e water gas production 
rate of both components was variable with time. iDello Russo et al.l ( 120071 ) observed a de- 
crease in the water production of B by a factor of 6 over a 6-day span of time, and a factor 
of 2 in a single day, 14.6-15.6 May 2006. The major "crumbling" event in B likely occurred 
near 1 April, closer in time to the Spitzer observations, which were likely even more affected 
by changes in water ice and vapor production. Their data also indicate variability in C, with 
a 50% increase (4.8cr) in a single day. In fact, although the water production rate was three 
times greater in B than C on 9 April, B's water production rate was half that of C's on 15 
April. Thus the relative total strength of the water bands in the two comets is a matter of 
the timing of the observations, and if the rate at which the emitting species remains in the 
extraction region of the spectrograph differs for the gas and the dust (which have different 
trajectories and speeds), the contrast of the water bands over the dust continuum will also 
simply be a matter of timing. 



4. Discussion 



Comets are composed of interstellar materials that have remained intact, material cre- 
ated during the formation and evolution of the protosolar nebula, and remnant material that 
was further processed in the nebula at a later time. Due to the presence of a variety of trans- 
port mechanisms, material that was once in the innermost regions of the solar nebula may 
eventually find their way into the comet-formation zone, where they may be incorporated 
into the comets themselves. The degree to which this occurred, its timing, and the nature 
of the processing will be imprinted on the comet population. 

The actual chemical content of the material in the early solar nebula is the result 
of a complex network of gas-gas chemical reactions, gas-grain reactions , dus t cond ensa- 
tion, annealing, and an d sublimation (see, for example the models of iGaill (120041 ) and 
Wehrstedt fc Gaill (120081 )). Because the temperature of the nebula will be highest close 
to the sun, the processing occurring there is likely to dominate, but imposes the condition 
that the radial transport must be efficient to introduce the material into the comet-forming 
region. 
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4.1. Processed Silicates 



The observed degree of crystallinity in SW3 and the other comets studied spectroscop- 
ically is well in excess of that predicted for the first million years by these models. This 
suggests that: the process by which crystalline material is produced is more efficient than 
these models predict, the radial mixing is more efficient than predicted, or it lasts much 
longer than 1 Myr, or other mechanisms for the production and transport of crystalline 
materials are also in play. 



Thermal processing of interstellar dust grains in the inner solar nebula (jChick fc Cassen 



19971 ) offers a natural source for the presence of crystalline silicates in the solar system. The 
resulting grain chemistry will be dependent on the availability of oxygen. Inside the water 



disso ciation line, OH will also react with C to make CO and CO2 (lGailll200ll 120021 : [Wooden 
20081 ). Low O fugacity favors the condensation of Mg-rich silicate crystals and Fe metal, 
and the annealing of amorphous Mg-Fe silicates into Mg-r i ch sil i cate crystals (forsterite - 
Mg2Si04, and enstatite - MgaSigOe) and Fe metal JGaillll998[ [20o3 : IWoodenll2008h . However, 
an influx of water from the outer disk can shift the balance toward the production of more 
Fe-rich crystalline silicates (fayalite - Fe2Si04, and ferrosilite - Fe2Si206), and less Fe metal. 

While the grain processing is occurring, radial mixing in the protosolar nebula can dis- 
tribute this processed material beyond 5 AU, into the reg ion where comets were assembled 
JCaill I2OO4I : iKeller fc Ga"iil l2004i IWehrstedt fc Gail! 120081 ). However, the balance between 
crystalline and amorphous material will have a strong radial dependence, as will the mineral 
content of the material. In the innermost edge of the surviving dust disk, crystalline forsterite 
dominates, but is quickly replaced by enstatite. Beyond 10 AU, these pure crystalline species 
make up only a few percen t of the silic ate material, which is now almost completely amor- 
phous. Early estimates by iGaill (J2004J ) of the resulting mineralogy near 20 AU placed the 
fractional composition of enstatite at ~13% and forsterite at ~1%. These estimates, however, 
were dependent on having complete chemical equilibrium between enstatite and forsterite, 
and mo re realistic modeling was expected to reduce the fraction of enstatite. In more recent 
models (IWehrstedt fc Gailll2008l ) these minerals have roughly equal abundances, but neither 
are present after 1 million years at more than a few % beyond 10 AU. 



Cieslal (120091 ) has modeled the transport of processed inner nebula dust by bipolar jets. 



where the dust decouples from the gas and is deposited onto the surfaces of the outer disk. 
The grains are them mixed into the disk by turbulent diffusion. The question remains 
whether the diffusion can mix the grains into the central plane, where comet growth is likely 
greatest, on a time scale short enough for that material to be included in the interiors of the 
comet nuclei. The tremendous degree of fragmentation of component B gives us confidence 
that we are seeing a substantial portion of what was once the interior of the comet parent 
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body, which probably accreted early in the life of the solar nebula, and crystalline material 
is certainly present, albeit not in the quantities seen in some other comets, such as HB. 

Simi l arly, r adiation pressure can also move material over substantial radial distances. 



Vinkovid ( l2009l ) has shown that this process can be quite effective if the radiation from the 
disk itself is also included. The disk radiations helps to levitate the grains, so that photons 
from the star cause the grains to glide over its surface to substantial distances. Whether 
this mechanism is capable of transporting processed material to the comet-forming region 
depends on the time scales of over which it operates effectively compared to the availability 
of processed material to be transported, and the time scale of comet formation. 



Harker fc Desch I (|2002| ) have suggested that nebular shocks at 5-10 AU can produce 
local temperature spikes that are capable of thermally annealing grains, obviating any need 
for large scale transport. However, their model predicts essentially no crystalline material 
should be found beyond 10 AU, and as such would be absent from comets formed in the 
trans-neptunian region, which is the origin of the ECs, such as SW3. 

Another possible mechanism for producing crystalline material closer to the comet- 
formation zone is throu gh heating by exoth ermic chemical reactions in organic refractory 



materials in the grains. iTanaka et al.l (|2010[ ) suggest that such reactions may be initiated 
at temperatures of only a few hundred Kelvins, and are sufficiently energetic to produce 
temperature spikes capable of annealing the grains. These would still need to be transported 
to regions in the solar nebula where the temperature has remained below 50 K. 

To summarize, a number of mechanisms exist that are capable of producing and trans- 
porting processed silicate material into the comet-forming zone. But it is unclear if any can 
produce the amounts of crystalline material seen in comets (even the relatively low crystalline 
fraction in SW3) actually observed. 



4.2. Elemental Abundances 

One of the aspects of solar system formation that we want to understand is the distri- 
bution of elemental abundances in various solar system materials, and how they are related 
to the bulk composition from which the solar system formed. In general, the "cosmic" 
abundances, which are really only those representative of stars formed at the galactocentric 
distance of the Sun, and at the present epoch, are derived primarily from spectroscopic anal- 
ysis of solar photospheric lines, and laboratory analysis of those meteorites thought to be 
most "pristine", least "processed", or "primitive" - the CI carbonaceous chondrites. It is to 
be expected that, as one considers materials condensed in regions of the solar system least 
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hospitable to the survival of the more volatile species, that those atoms significantly locked 
up in volatile materials wo uld be dep l eted compared to the Sun, having been vaporized, and 
indeed many of them are. lAsplundl (120091 ) summarizes the state of our current knowledge 
of the abundances of both the Sun and the CI chondrites. Of particular significance to our 
discussion here are the abundances of C and O, which are major atomic species in potentially 
volatile materials, and Si, the element to whi ch meteoritic abundances and some cometary 
abundances are normalized to (lAsplundll200^ 



Over the years, abundance determinations have evolved, particularly the photospheric 
ones, due to the difficulty of deriving true abundances from spectroscopic line analysis. For 
Si, the abundance is currently known to ±0.04 dex, and has not changed by more than 
~0.02 dex over the past decade or so. By contrast, the abundanc es of C and O have changed 
significantly since the classic work of lAnders &: Grevessd ( 119891 ). Much of the revision was 
driven by a problem dating back decades, when the depletion of elements from the interstellar 
gas was being addressed. 

The first major analyses of the atomic abundances in diffuse interstellar clouds that in- 
cluded species most easily studied with ultraviolet transitions that became available with the 
launch of the Copernicus observatory. The result was that mo s t spec ies se emed to be depleted 
relative to solar values (see the reviews of ISpitzer &: Jenkins! (Il975l ) and lSavage fc Sembach 
(119961 )). In some cases, such as Al and Ca, the depletions in the gas phase with respect 
to solar were a factor of ~10^, while in others, such as C and O, it was between 3 and 5 
times. As the measurements became more refined and more lines of sight were sampled, 
this pattern remained generally the same, with slightly different depletions along different 
lines of sight. The depletion seemed to be correlated with the condensation temperatures of 
the main sinks of the elements under circumstellar conditions originally, but later analyses 
were done within the framework of de nse outfiows f r om ev olved stars, from which the grains 
originally condensed (see for example ISavage et al.l (119921 )). Throughout, these abundances 
were compared to the most current solar values. Despite the relatively small depletions of C 
and O in these diffuse clouds, they were important because, due to their overall large cosmic 
abundances and importance in grain composition, any deficit or excess in the depletions 
could impact the availability of the elements n eeded to make the c orrect amount of dust to 
reproduce the observed interstellar extinction (ICardelli et al.lll996l ). 



^Results of the impact ionization mass spec t rometer data on IP /Halley have b een presented in terms 
of normahzations to Si (jLangevin et al.l Il987l : iLawler et al.l Il989l ) and to Mg (jjessberger et al.1 Il988t 
Jessberger fc Kissell 1991). The use of isotopic ratios to identify ion abundances was stated to be bet- 
ter for Mg by iJessberger et al.l (|l988f ) , so we have plotted both. See the Appendix for a more complete 
discussion. 



-21 - 



At the same time it was realized that t he photospheric abundances derived for the Sun, 
such as those of lAnders fc Grevessd (119891 ) and earher works, were inconsistent with those 
derived from youn g stars being born toda y, with the Sun being overabundant in C and O 
by a factor of ~2 (jSofia et aLlll994l . Il997l ). If the "cosmic" abundances were smaller, the 
discrepancy would be minimized, but would also lower the amount of C and O available to 
make dust. 



Since the work of lAnders fc Grevessd (119891 ) the solar photospheric abundances of C and 
O have declined significantly. For C, th e log of the abu ndance (whe re H is defined as 12. 00) 



is, depending on the source, 8. 43 ±0.05 (lAsplundll2009l ). 8.39±0.04 (iLodders et al 



8.78±0.07 (ICaffauet al.ll2010f ). compared to 8.9 3±0.0435 lAnders fc Grevess^ 



2009a 



2009a 



2009a 



b) 



or 8.5 0±0.06 (ICaffauet al.ll2010f ). down frorn the value of 8.56 ±0.04 of 
( Il989h . F or O, it is 8 69±0 05, JAsplundl l2009[ ). 8.73± 0.04 Jhodders et al 



Anders &: Grevesse 



2009a 



b|), or 



1989). Thus 



the ne wer values have dropped b y about 13% ( Caffau et al.l 120101) to 32% (ILodders et al. 



2), and O has dropped 33% JCaffau et al.lboioh to 47-48% JAsplur^bood : Ihodders et al 



^. 



It is with these revised "cosmic" abundances that modern comet and meteoritic results 



must be compared. For simplicity, we will use the values from lAsplundl (120091 ) for our analysis 
of SW3. 



4-2.1. Caveat Emptor 

Two issues that are easy to ignore in the spectroscopic analyses of comet abundances are 
the possible presence of major sinks of elements in undetected species, and the unknowing 
removal of species from the sample. The first can lead to a "depletion" that does not exist. 
The second to assign a real depletion to the wrong cause. 



TielensI ( 120051 ) notes that the abundance of C and O drop further from their values in 



the diffuse interstellar clouds discussed previously to denser molecular clouds - the sites of 
star formation, and the starting point of the solar nebula. It is possible that much of the 
O becomes locked up in molecular O2 which is not easily measured, while the source of the 
increased deficit for C remains unclear. The increasing depletion of C and O is independent 
of any solar abundances, and serves as a warning that both C and O may be stored in forms 
that are not immediately detectable in spectroscopic studies. Species so stored may or may 
not later become processed into detectable forms. 

The C and O may also be preferentially removed from the comet-forming regions long 
before the comets are assemble. CO (and CO2) acts as a sink for both species. As dust 
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settling to the mid-plane of the nebular disk proceeds, it leaves behind uncondensed gases, 
such as CO, that will not be incorporated into the comets themselves. Recent observations 
of young protostellar disks with ages of a few Myr are consistent w ith the gas remaining in a 
thick flared disk after the dust has largely settled ( lAcke et al.ll2010l ). Bodies being assembled 
there will be deficient in those elements locked into the major gas species - H2 and CO. In 
this case, it is not simply a matter of the volatiles being too warm to condense, but being 
significantly separated spatially from the refractories. There is therefore no a priori reason 
to expect the C and O abundances in comets to be exceptionally close to solar to begin with. 
But comparing them to the abundances of other solar system materials has the potential of 
shedding light on the evolution of the solar nebula. 



4-2.2. Rock-Forming Elements 



Figs. US] and [T3] show the abundances of a number of elements relative to both Si 
and Mg in SW3 and three other comets studied spectroscopically with either Spitzer or the 
Infrared Sp ace Obse r vatory (ISO), normalized to the solar photospheric abundance ratios, 
taken from lAsplundl ( 120091 ). For the majority of the rock- forming elements, SW3 is similar 
to the other comets, and close to or somewhat less than solar. For comparison, we also 
show the abundances for IP/Halley, based on the impact ionizat ion mass spectro r neters 
of t he Giotto and VEGA space craft based on two studie s , that of IJessberger et al.l (Il988[ ) 
and I Jessberger fc Kisselll ( 1l99ll ). and that of lLawler et al.l (19891). Als o shown are the mean 
ab undances for th e CI carbonaceous chondrites, from iLodders et al.l ( l2009al Jbl) (also listed 



m 



Asplundl ( I2OO9I )). For Mg (in the Si- normalized plot). Si (Mg- normalized), and Fe (both 



normalizations), the abundances of all five comets studied spectroscopically, including both 
components of SW3, are consistent with both the CI and solar photospheric values. For Al, 
S, and Ca there is considerable scatter around the the CI and photospheric values for SW3, 
but the uncertainties are larger, as these elements are present only in minerals with trace 
amounts (diopside, smectite nontronite, and niningerite) in the models. The abundances of 
these elements in these five comets are also consisten t with the m a ss spe ctrometer data on 
Halley, although the uncertainties for Ca are large in iLawler et al.l (119891 ). and they did not 
report an abundance for Al. 



4.2.3. The Oxygen Abundance 



For all of the objects shown - the five comets studied spectroscopically, both sample 
results on Halley, and the CI chondrites, O is depleted with respect to the solar values. The 
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spectral studies specifically include H2O vapor and ice as components, while the Halley data 
do not, and the chondrites are not as water-rich as the comets. The depletion of O in the 
comets determined from the spectral models might result from: 

• the water model underestimating the actual amount of H2O vapor present 

• the spectral model underestimating the H2O ice abundance 

• the dynamics of the gas compared to the entrained dust causing the H2O to be 
underestimated compared to the dust (Si, Mg, etc.) 

• O being sequestered in some dust species to which the spectral models are insensitive 
(large or spectrally featureless grain materials) 

• O being in gaseous species not detected in the spectra 

• O being truly depleted in comets 

The H2O vapor abundances are determined by the bands in the 6 /zm region, which are 
well-defined, and would need to be off by a factor of ~3 (regardless of which normalization 
is used) to bring them into agreement with solar values. 

The H2O ice band occurs at the minimum between the 10 /xm and 20 /im silicate bands 
and is not easily masked or blended with other features. Very large ice particles that were 
missed in the particle size distribution would be spectrally weak and might go undetected. 
The likely presence of large chunks of ice in component B might explain its lower spectrally- 
detected abundance than in C. But C is also under-abundant in O, so this cannot be the 
sole answer either. 

For a fixed spectral extraction size, faster-moving species, such as gas molecules, will 
travel out of the extraction zone than the more slowly-moving dust grains. For a grain 
ejection speed imparted by the outflowing gas, 

Vej ~ 2.3{pacmy^^'^ms'^ (3) 



(iReach et al.ll2010l ) a grain with radius a ~ 1 /im and p ~ 2 g cm^ will be moving with v ~ 
0.3 km s^^, or about 1/3 the speed of the gas. So compared to the dust grains responsible for 
the spectral features modeled, the water gas will be under-represented, which would cause 
O to be under-represented. Note that this also would suggest that larger spectrally gray 
grains may be over-represented in the extraction aperture, as they lag the smaller grains. 
This may influence the PSD derived from the spectral fitting, which is further complicated 
by the presence of multiple "nuclei" in the extraction aperture for fragment B. 
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The Stardust sample of comet Wild 2 contains dust material that does not fall into 
the minerals used in the spectral fitting. The largest grain first examine d in detail was 



one s imilar to calcium-aluminum-rich inclusions (CAIs) found in meteorites (jZolensky et al. 



20061 ) . which contain many oxygen-bearing minerals. The Stardust grain contains spinel 
(MgAl204), anorthite (CaAl2Si208 in its pure form), and gehlenite (Ca2 A^AlSiO?] ) . Othe r 
meteoritic CAIs are generally rich in Si02, AI2O3, MgO, and CaO groups (JKrot et al.ll200ll ). 
Such mineral species, especially if they were in 1 mm sized grains, would be missed in 
the spectral modeling. A comparison of the Gi otto and VEGA probes showed that the 
scattered light was dominated by mm-sized grains (JMcDonnelllll99ll ) . In SW3 itself, analysis 
of the Spitzer Miltiband Infrared Photometer (MIPS) images of the dus t tails and dust trails 
indic ates that the "bump" seen in the PSD continues to cm-sized grains (IVaubaillon &: Reach 
2OIOI ). This is supported by the radar measurements (JHowell et al.ll2007l ). So the contribution 
by large grains to the O abundance remains uncertain, but may be significant. 

In principle, CO and CO2 gas might be another possible sink for the O not detected in 
the spectral mod ehng. CO is known to be severely depleted in SW3 in particular, compared 
to other comets (JDiSanti et al.ll2007l ). If CO is a major source of 0-depletion from the dust, 
it is not in the form of cometary CO in SW3. There is tremendous scatter in the abundance 
of CO amongst the comet population (see below) and much of this may be the result of the 
separation of the dust and gas in the solar nebula, similar to that seen in young protostellar 
disks, as discussed previously. 

Another source of the depletion of O in the spectral studies is its incorporation into 
organic species. In the Halley samples, the bulk of the C is in for form of complex organics 
- the "CHON" mate rial JPomenkova et "aDll992l . Eooi : lFomenkovalll999h . For the 40 CHON 
particles studied by IJessberger et al.l (119881 ) , the rati o of O to C atoms for the sample has 
considerable scatter ranging from 0.2 to 10. In fact, IJessberger fc Kisselll (1199 ll ) concluded 
that "most of the oxygen is not from the silicates, but from an organic phase." Doubling the 
O abundance in the spectrally analyzed comets would bring their values of O closer to that 
of the CI chondrites, but probably not to solar values. The O stored in organics is included 
in the Halley sampling, but not the O stored in H2O, as the sur vival of ice particles at the 
distance of the Halley sampling is thermodynamically excluded (IJessberger et al.lll988l ). A 
gas/dust ratio of ~2 would bring the Jessberger et al. abundance close to solar, while the 
ratio would have to be considerably larger, ~5 or more, for the abundance determination of 
Lawler et al. Thus both the in situ sampling and spectral models of comets may be missing 
significant 0-bearing material, but a different material in each case. 

In summary: 

• O seems to be depleted in all comets, including SW3, compared to solar values. 
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• some of the under-abundance may be due to under-sampling the O tied in H2O gas 
in the spectral models, which exits the sampling beam faster than the dust. 

• some of the under-abundance in the spectral models may be due to not including 
0-rich material that is tied up in large rocky grains and CHON material. 

• the in situ measurements of Halley, which include some large grains and 0-bearing 
CHON, are still somewhat under-abundant in O, but the O in the form of H2O is missed. If 
the dust-to-gas ratio is one or more, the results approach the solar values. 

• A lot of O was likely tied up in CO gas in the nebula which remained "at altitude" 
while the dust settled into the mid-plane of the solar nebula. Thus most of the objects in 
the solar system may have formed under conditions where H, C, and O were depleted. 



4.2.4- The Carbon Abundance 

One of the most striking aspects of the model fits to the Spitzer spectra of SW3 is 
the high abundance of amorphous carbon. While the C abundance in most comets is sub- 
solar, and similar to that of carbonaceous chondrites, it appears near to solar in SW3, 
within the uncertainties, which are somewhat larger t han the othe r objec ts studied. It is 
also considerably high er than the Hall e y ana lysis by iLawler et al.l ( 1l989l ) but marginally 
consistent with that of iJessberger et al.l ( 1l988l ) which is at least a factor of two higher than 
that of Lawler et al. 

Using an independen t set of observations of SW3 with the Michelle spectrograph on the 
Gemini Nort h teles cope, iHarker et al.l (I2OIOI ) modeled both components of SW3 with the 
Harker et al.l (120021 ) thermal emission code, and found that the the relative abundance of 
amorphous C in the anti-sunward coma was between ~ 45-23 wt% for B and ~60-42 wt% for 
C. Using a similar treatment, the amorphous C abundance was determined to be 21 wt% for 
HB, 15 wt% for C/2001 Q4 (NEAT), and 28% for Tl. The fractional abundance of carbon 
compared to silicates decreased with increasing distance from the nuclei of both components. 
The median values within the ranges found from the Michelle data are significantly higher 
for SW3-B that either HB or Tl, while SW3-C was even higher. Thus these independent 
analyses, using different mineral input spectra, a different modeling code, and observing the 
SW3 components at a different time than the Spitzer observations, yielded similar results 
for the high abundance of C in SW3 compared to other comets. 

At th e same time, SW3 is highly deficient in alraost all C-bearing raolecules except HCN 



and CO2 fIVillaneuva et al.ll2006l : lOiSanti et al.ll2007l : lOello Russo et al.ll2007l : iLis et al.ll2008 
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Reach et al.ll2009l ). Carbon combustion in the inner solar nebula is expected to result in the 
removal of amorphous C and its conversion to CO, CO2 and various C-bearing hydrocarbon 
molecules such as CH4 and C^H^, although ot hers such as C2He and CH3OH r equire a differ- 



ent so urce (lGailll2002 : 



Wehrstedt fc Gailll2008[ ). CO was detected in SW3-C bv lOiSanti et al. 



(I2OO7I ) with an abundance of C O/H2O = 0.5% ± 0.13% , making it at the lowest end of of 



the observed ratio in comets. iDello Russo et al.l ( 120071 ) reported that SW3 was depleted 
in most carbon-chain molecules, compared to the majority of other comets studied so far. 
The high abundance of amorphous C and simultaneous low abundances of CO and carbon- 
chain organics are all consistent with SW3 being formed in a region that did not inherit 
the products of C-combustion in the inner solar nebula. Further depletion of gaseous C- 
molecules compared to solid amorphous C will occur as dust sedimentation coours in the 
disk mid-pla i ie, as was discussed for previously for atomic O. The one discrepancy is CO2. 
Reach et al.l (120091 ) . using Spitzer imaging data, derived CO2/H2O ratios of 10% and 5% for 
B and C, respectively, which indicates that CO2 was not depleted. 

A significant (if not well-quantified) fraction of cometary CO may come from a "dis- 
tributed" coma source that might include refractory material, or f rom a complex network 
of chemical (and photo-chemical) reactions (jPierce fc A'Hearrul2010l ). and might behave dif- 
fere ntly. The CO abundance varies across the comet population by as much as a factor of 
40 (JMumma et al.ll2003t iBockelee-Morvan et al.l 120041 ) and as yet does not seem to be well- 
correlated with other species. A s i gnific ant complication arises from the CO band optical 



thickness. iBockelee-Morvan et al.l ( 120101 ) have recently analyzed both the CO J(l-O) and 



J(2-l) lines in Hale-Bopp, and have concluded that the CO f =1-0 rovibrational lines often 
used to detect CO - and provide evidence for a distributed source - are severely optically 
thick for heliocentric distances of 1.5 AU or less, and that in fact there is no evidence for a 
distributed CO source in this comet. 

One added uncertainty in the determination of the elemental carbon abundance is our 
uncertain knowledge of the amount of metallic Fe present in the dust, in the form of bulk Fe 
or FeNi, which could contribute to a featureless continuum in the models, and be responsible 
for some of the continuum attributed to amorphous C. Metallic Fe is a well-known component 
of chondritic porous interplaneta ry dust particles (IDPs), such as t he GEMS (Glass with 
Embedded Metals and Sulfides) JBradlev et al.lll992l : [Bradlevlll994h . ^ 



Finally, scattered solar light will contribute most strongly at the shortest wavelengths 



lOWhile the GEMS themselves possess infrared spectral features resembling the amorphous dust of 
the interstellar med ium, the IDPs in which they are embedded are similar to those observed in comets 
(JBradlev et al.lll999l) . The presence or absence of GEMS material in the Stardust samples is complicated by 
the melting of the aerogel in the high-speed capture process. 
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precisely where the models are most sensitive to the carbon content. We have placed an upper 
limit on this contribution at 11% at 5 ^m. (Appendix). This would reduce the amorphous 
carbon contribution by at most 10%, not substantially affecting our compositional results. 

In summary, the abundance of amorphous carbon derived from the spectral models 

• appears equal for B and C, within the uncertainties 

• is higher than the majority of other comets studied spectroscopically (in both inde- 
pendent studies) 

• is consistent with it a low abundance of CO and carbon-chain organic molecules if 
SW3 did not incorporate significant amounts of C-combustion products from the inner solar 
nebula 

In terms of the atomic abundance (not the form that it is in): 

• carbon appears nearly solar in both components of SW3, although the uncertainties 
are large. 

• carbon is somewhat higher in abundance than other comets whose composition is 
determined from spectral decomposition, and possibly higher than the in situ measurements 
of Halley. 

• the slightly higher atomic C may result from it being sequestered in a moderately 
refractory solid phase, instead of highly volatile ices 



4.3. Comet Homogeneity 

If the mineral content of the region where cometesimals grew and were assembled into 
their final cometary form occurred in a region where the composition of the dust mate- 
rial was evolving with time due to the delivery of processed material from the inner disk, 
we would expect that individual cometary nuclei would be chemically and mineralogically 
heterogeneous. Their surfaces would then have a different mineralogy than their interiors. 

The fragmentation of SW3 exposed material that was once deeply buried in its interior, 
and the further disintegration of component B assured a greater exposure of this material 
than in C alone. Both components have approximately the same carbon content (based on 
independent observations and spectral models) and degree of silicate crystallinity, suggesting 
that they formed from the same mix of materials. If B is overall exposing more deeply-stored 
materials than C, this would suggest that the components were assembled at a time when 
"salting" of the refractories by materials processed in the inner solar nebula was minimal. 
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There is no "time signature" evident. 

Given this argument and the compositional similarity of fragments B and C in the more 
refractory materials, as well as gaseous species (JDello Russo et al.l 120071 ). we may conclude 
that there was a substantial homogeneity in the pre-brea kup nucleus of SW3 . Furthermore, 
if we consider that the radius of fragment B is <0 .4 km (JHowell et al.l 120071 1. and fragment 
C is ~0.7 km (JToth et al.l 120051 : iHowell et al.l 120071 ). then we may consider that perhaps the 
entire pre-breakup nucleus is homogeneous throughout. 



5. Conclusions 



The inner solar nebula possessed the thermal and chemical ingredients for the production 
of silicate material significantly more crystalline than the raw interstellar dust from which 
it formed. Gas-gas and gas-grain chemical reactions likely converted some portion of the 
primordial amorphous carbon in ISM grains into a variety of carbon-containing molecules. 
Both the crystalline material and the amorphous C-depleted and hydrocarbon molecule-rich 
material were mixed radially, although cometary material formed very early or far from the 
inner nebula would not be affected to the degree of those formed later or closer to the sun. 
Spectroscopic observations of both B and C components of 73P/Schwassmann-Wachmann 
3 reveal a low degree of grain crystallinity compare to other comets (but higher than the 
ISM), higher amorphous carbon content than most other comets, and depletion of most C- 
rich molecules relative to other comets. These are all indicative of being more primitive than 
other comets studied remotely, and suggest it was formed earlier than most, or further from 
the warm inner nebula, or both. 

No significant difference in the crystallinity nor amorphous carbon content was observed 
between the two components in either the Spitzer data nor the independently-observed and 
modeled Michelle spectra by iHarker et al.l (l2010l ). The difference in the rates of nucleus 
fragmentation suggests that we are likely digging deeper into the pre-fragmentation nucleus 
in B than in C. Despite this, they appear spectroscopically very similar, arguing against 
significant heterogeneity and/or layering. 

One of the most significant differences observed, the relative strength of the water vapor 
bands, may not point to any intrinsic difference in the two components, but instead to their 
being significant temporal changes in their water vapor production, or significant differences 
in the distribution of water ice sources in the extraction beam. 

All comets studied spectroscopically exhibit a general under-abundance of 0-rich species. 
It is likely that much of the O is tied up in grains that are spectroscopically featureless, or 
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at least not easily detectable with current techniques. These do not, however, explain the 
deficit in the O abundance when measured in situ (such as Halley), which is likely dominated 
by the exclusion of the O on the form of H2O gas. Some of it may be in the disk gas that 
that remains at high altitude in the proto-solar nebula as the dust from which comets form 
settle to the mid-plane, as has been seen in other protostellar disks. 

A better understanding of the chemical and physical evolution of the early solar system 
through the analysis of cometary material will require many studies using a variety of com- 
plementary techniques, both remote sensing (spectra, imaging) and in situ sampling (dust 
fluence experiments, sample returns, remote sample analysis). 
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A. Ground-Based Observations of 73P/Schwassmann-Wachmann 3 

In addition to observations of 73P/Schwassmann-Wachmann 3 obtained with the Spitzer 
Space Telescope we also carried out ground-based observations of components B and G 
using the Michelle spectrograph on the Gemini North telescope and on NASA's Infrared 
Telescope Facility (I RTF) using the Spe X and BASS instruments. The Michelle results are 



discussed elsewhere (JHarker et al.ll20ld ). Here we use the IRTF observations to provide a 
second estimate of the degree to which scattered solar radiation contaminated the Spitzer 
spectra, determine an estimate of the albedo of the grains, and derive the dust flux parameter 
A{6)fp commonly used in comet studies. The timing of the Spitzer, BASS, and Michelle 
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observations, relative to the light curves of the comet fragment^lj, is shown in Fig. [151 



A.l. SpeX 

SW3-C was observed on 20 April 2006 (UT) using the SpeX spectrograph on NASA's 
Infrared Telescope Facility. Observations were made using the cross-dispersed echelle grat- 
ings in both short- waveleng th mode (SXD) cov ering 0.8-2.4 yum and long-wavelength mode 
(LXD) covering 2.3-5.4 fan (JRayner et al.ll2003l ). All observations were obtained using a 0.8 
arcsec wide slit, corrected for telluric extinction and flux calibrated. Because no suitable 
AOV star was observed close enough in airmass and hour angle to the comet observations, 
we used data obtained on HD 25152 on a later date with similar sky c onditions, airmass 



and hour angl e. The data were reduced using the Spextool software ( IVacca et al.l 12003 



Gushing et al.l l2004j ) running under IDL. Details of the observations are listed in Table 4, 
and the resultant spectrum is shown in Fig. Uni In this figure, the dropping flux density at 
the shorter wavelengths, followed by the rise at longer wavelengths, illustrates the transition 
from scattered solar radiation to thermal re-emission by the cometary dust. 



A. 2. BASS 

Both B and C components were observed with The Aerospace Corporation's Broad-band 
Array Spectrograph System (BASS) on 18 and 19 May 2006 (UT). BASS uses a cold beam- 
splitter to separate the light into two separate wavelength regimes. The short-wavelength 
beam includes light from 2.9-6 fim, while the long-wavelength beam covers 6-13.5 fim. Each 
beam is dispersed onto a 58-element Blocked Impurity Band (BIB) linear array, thus allow- 
ing for simultaneous coverage of the spectrum from 2.9-13.5 fira. The spectral resolution 
R = A/AA is wave length-dependent, ra nging from about 30 to 125 over each of the two 
wavelength regions (JHackwell et al.lll990l ). The details of the observations are listed in Table 
5 and the spectra of components B and C are shown in Figs. [T7| & [TSl 

The BASS spectra of B and C are nearly indistinguishable, with the exception that the 
relative strength of the silicate band is slightly stronger in C than in B (Fig. [T^ . This is 
consistent with what we observed in the Spitzer IRS data. 



^^ Adapted from Sechii Yoshida's comet web site |http://www. aerith.net/index.htmll and used with permis- 



sion 
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A.3. BASS CCD 



On 19 May 2006 (UT) we also obtained observations of SW3-C using the CCD guide 
camera of BA SS and the Blue Conti nuum (0.44 /xm) and Red Continuum (0.71 /xm) "Hale- 



Bopp" filters (iFarnham et al.l ( 120001 ): hereafter Be and Re)- Observations were calibrated 
using the flux calibration stars (HD 186408, HD 164852, and HD 170783) for the filters. 
Observations of 73P-B were unsuccessful due to bright twilight. The details of the CCD 
observations are listed in Table 6. In addition to the comet filters, we obtained a deeper 
image using a broad-band Cousins-Kron R filter, but this was not used to extract photometric 
information. Photometric information was extracted from the Be and Re images a circular 
software aperture matching the entrance aperture of the BASS dewar. 



A. 4. Scattered Light Contribution to the mid-IR Spectra 

In Fig.[20]we show the 3-13 /im spectrum of component C obtained on 10 May 2006 (UT), 
along with the Bq and Re flux densities obtained the same night. To this we added the SpeX 
data, scaled to roughly match both the CCD photometry and mid-IR BASS spectrum. The 
match in the 3-5 /im region where the SpeX and BASS spectra overlap is only approximate, 
as they were obtained when the comet was a different heliocentric distances (SpeX - 1.16; 
BASS - 0.98), and hence slightly different temperatures. 

Also shown is a slightly reddened smoothed solar spectrum that approximately fits the 
combined CCD and SpeX data, plus a simple grain emission model. The model consists 
of mixture of four grain components: amorphous olivine and pvroxeny^l. crvs t alline olivine. 



and amorphous carbon. Optical constants were taken from iDorschner et al.l (119951) for the 



amorp hous silicates with Mg/Fe = 1.0, that of the crystalline forsterite from iKoike et al. 



( I2OO3I ). and amorphous car bon (specifica l ly arc discharge from amorphous carbon electrodes 



in an Ar atmosphere) from IZubko et al.l (119961 ). A Mie scattering model of the three amor 



phous constituents was calculated for solid spheres using size dis tribution for IP/Halley 



determined with the DIDSY instrument on the Giotto spacecraft (iMcDonnell et al.lll987| ) 
with size cutoffs at 0.06 and 70.0 /im. To this was added the crystalline forsterite, and the 
temperature (for this purpose, assumed to be the same for all the components) adjusted 
to fit the BASS data. The purpose was not to obtain an independent assessment of the 
grain composition, but to approximate the scattered light and thermal emission in order to 



^^Strictly speaking, the terms "olivine" and "pyroxene" refer to pure crystalline end-grroups of the min- 
erals, but we use them here instead of "amorphous silicates of roughly olivine and pyroxene stoichiometry" 
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obtain two quantities: the grain albedo, and the degree that scattered hght contributes to 
the wavelengths included in the Spitzer spectral modeling. Here we use the fact that B and 
C are spectroscopically nearly identical, so that the merged data and model for C can be 
used to estimate the scattered light contamination for both B and C on the dates of the 
Spitzer observations. 

For component B a simple extrapolation of the SpeX data short-wavelength spectral 
slope to longer wavelengths would suggest 1% for the scattered light at 5 /im, but looking 
globally at the SpeX+BASS observations would imply that it is closer to 3%. For C the 
fraction of the light at 5 /im that is due to scattered solar photons is probably about 7% for 
the da ta obtained on Mar ch 17 (when the dust was cooler than that of Spitzer observations 



of B). iReach et al.l ( 120091 ) put it at 6% in the 4.5 /xm IRAC image, obtained on March 28. 



While the solar spectrum used is slightly reddened to approximate the flux levels in the 
Hale-Bopp filters, the value of 7% assumes that t he gr ains are becoming a little bit "blue" 
in the infrared, consistent with IJewett fc MeechI (119861 ). If the grains are "gray" scatterers 
at all wavelengths, the number is higher, and possibly closer to 10.7%. The scattered light 
component will most heavily affect the grain component that dominates at the shortest 
wavelength, which in our models is the amorphous carbon. However, even a 10% change in 
the derived abundances does httle to change the fact that the carbon abundance in the dust 
of SW3 is higher than the other comets to which the same spectral models has been applied. 



A. 5. Grain Albedo 



For this simple model, the bolometric albedo (net integrated scattered light divided 
by the sum of the scattered light and thermal emission) was ~0.13. At the phase angles 
that these observations were made (60°-80°) this is somewhat low er than many coraets, a l- 
though similar to the curve for C/1980 Yl (Bradfield) from Fig 1 of iKolokoloya et al.l ( 120041 ). 
Note that this definition of the albedo is the same as that o f iGehrz fc NevI (119921). s ome- 



times called the bolom etric albedo. This differs from that of iHanner fc NewburnI ( 119891 ) and 
Tokunaga et al.l ( 1l986l ). which is the geometric albedo, usually multiplied by a normalized 
phase function. The latter are generally mu ch smaller than the bolometric albedo. For a 
more complete discussion of these terms, see iHanner &: NewburnI (119891 ). 
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A.6. Variability 

We observed both components with BASS over two nights. The observations of B and 
C were interleaved with each other, with observation of flux calibration stars, and with the 
CCD imaging. These provided three sets of data on B spanning 5 hours on 18 May, and a 
single set of observations on 19 May. For C we had two such sets on 18 May spanning 4.5 
hours and two sets on 19 May spanning 1.7 hours. 

For the 18 May observations of B, the resultant flux levels were identical to within the 
uncertainties, despite the fact that the second set of observations were obtained after the 
start of twilight, and the third many hours after sunrise. In both cases, the centering of 
the comet in the BASS entrance aperture was frequently checked by adjusting the centering 
to maximize the signal. The simplest explanation for the agreement of the three data sets 
is that B was stable during these observations. It is unlikely that errors in centering and 
tracking were precisely cancelled out by source variability, and gives us confidence that the 
observed changes in brightness in our daytime observations of either comet are likely real. 
The 30% drop on brightness on the following night seems large at first sig ht, but such rapid 



changes were observed in the H2O production rates in B reported by iDello Russo et al. 



(120071 ). Our BASS observations were obtain ed closer to the major breakup event of B than 



the observations of iDello Russo et al.l (120071 ) . when large variations might be expected. 



Contrary to our expectations, C seem to exhibit greater changes in mid-IR emission than 
B. If a significant fraction of C were "old" surface that survived from before the original 1995 
breakup, rotational modulation of "old" and "new" surface to solar heating might explain 
this behavior. A similar behavior is seen in the data obtained at shorter wavelengths. 

Optical observations of comet dust can be expressed in terms of the product A{9)fp, 
where A{9) is the grain albedo as a function of phase ang le, / is the dust fill ing factor. 



and p is the radius of the aperture used to measure the flux (jA'Hearn et al.lll984l ). typically 
expressed in units of cm. For comae in an isotropic steady-state outflow, or a coma that 
would otherwise follow a p~^ surface brightness profile, A{6)fp is independent of p. A{6)fp 
is used as a proxy for comet dust production rates, but also depends on the comet's grain 
size distribution, and dust scattering properties. 

For SW3-C, the flux density in the Re filter was Fa=4.28 x 10"^=^ W mr'^pmr^ (AFa=3.05 
X 10~^^W m~^; Yy=72.b mJy), and A{6)fp= 237 cm. The aper tures used were V!7 (123 km) 



in radius, equivalent to the BASS mid-IR entrance aperture. iBertini et al.l (120091 ) present 
Johnson i?-band A{6)fp values measured in a 4000 km radius aperture throughout April and 
May 2006. Our Rq filter image does not have the signal-to-noise ratio required to properly 
measure A{6)fp in a 4000 km aperture. To properly compare our Re value to their values. 
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we measured the azimuthally averaged surface brightness profile in our broad-band J?-band 
image to determine a scale factor (0.51) for ou r p = 123 km measu rement. We compare our 
corrected value of A{6)fp=122 cm to those of iBertini et al.l (120091 ) in Fig. [511 



To verify fragment C's variability, we examined the A{6)fp data in Fig. [21] A sub- 
stantial day-to-day variability is app arent in the A(9)fp d ata, although i t dep ends on the 
interpretation of the errors given in iBertini et al.l (120091 ). iBertini et al.l (120091 ) state that 
their A{6)fp errors are on the order of ±10% (we have used 10% in Fig. [2ID, and that these 
are primarily flux calibration errors. If the flux calibration error is independent of the night 
the comet was observed, then the relative scatter in A{6)fp may be real, and consistent 
changes seen in the BASS data. 

Such modulation would not normally be apparent in a comet coma because short 
timescale variability would be averaged out by large physical apertures. However, comet 
SW3 passed very close to the Earth, and was 0.09~AU away at the time of the BASS ob- 
servations. Photometry in a 1'.'7 radius aperture is sensitive to 123f ~^ s timescales, where v 
is the outflow speed in km s~^ of the dust being measured. Micron and smaller sized dus t 
grains are ejected at speeds of order 0.1 km s~^ at 1 AU from the Sun (JLisse et al.lll998l ). 
which corresponds to 20 minute timescales. Thus, dust comae variability on rotational pe- 
riod timescales could be easily observed, especially given the small grain dominated grain 
size distributions of comet 73P (dn/da oc a~^'° and a~^'^, §3.1.3). 



B. Comparison to Halley in situ Results 

B.l. Sampling of IP/Halley 

In order to provide a comparison to comet abundance data obtained through direct sam- 
pling rather than spectral modeling we have used the results of the three impact-ionization 
mass spectrometers flown through the coma of IP/Halley on the Giotto spacecraft (PIA 
instrument) and two VEGA spacecraft (PUMA-1 and PUMA-2 i nstruments) . T he main 



results of abundances deriv ed from these instru ments are those of iLangevin et al.l (Il987l ) , 
Jessberger et al.l (119881 ). and lLawler et al.l (119891 ). Due to the limited data transfer rates, the 
bulk of the spectra were transmitted in compressed modes, while one out of every 30 were un- 
compressed. In both cases, ions with energies up to either 150 eV or 50 eV were measured , 



and these comprised the short and long spectra, respectively (I Jessberger fc Kisselllll99ll ). 



op 

lL 



Langevin et al.l (Il987l ) inclu ded 2284 compressed spectra, combined from all three instru- 



ments. iLawler et al.l (119891 ) began with the same spectra but removed those sets plagued 
by excessive noise, high backgrounds, and spectra with double peaks. They carefully ex- 
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amined the data sets for solar i sotopic ratios, which helped to eliminate contaminating 



effects, a procedure not used by iLangevin et al.l (Il987l ). The result was 433 compressed 



PUMA-1 spectra (284 short and 149 long), 29 uncompressed PUMA-1 spectra, and 32 un- 
compresse d PIA spectra. Howev er, they only tabulated the results of the 433 compressed 
data sets. IJessberger et al.l (119881 ) . used the 74 uncompressed PUMA-1 spectra. These were 
also checked for reasonable isotopic ratios. W hile these last data may be the cleanest, they 
represent a sample size only one-sixth that of iLawler et al.l (jl989[ ). 



Jessberger fc Kisselll (Il99ll ) took the uncompressed spectra of the ions and made an ap- 



proximate conversion into total atoms by estimating the total yields of th e instruments. How- 
ever, t he conversion was probably only accurate to within a factor of two. IJessberger &: Kissell 
(|199ll) tabulated th e ion ratios of Si, S, and Fe to Mg their PUMA-1 data with that of 
Lawler et al.l ( 119891 ). and the ratios for the ions are quite similar, even to the relative ratios 
in the short and long spectra separately. 



B.2. Normalization to Si versus Mg 

Throughout the literature on cosmochemistry, various groups have determined the rel- 
ative abundances o f the elements by norm alizi ng their values to (^usually) either Mg or Si. 
The abundances of IJessberger et al.l (119881 ) and IJessberger fc Kisselll ( Il99ll ) used Mg, as the 
isoto pe ratio discr i minan t used to select the cleanest data were considered to be better for 
Mg. iLawler et al.l (Il989l) used S i for t heir own ion data, and converted the yield-corrected 



atomic data o f I Jessberger et al.l (1l988l ) into Si-normalized values. The CI chondrite data of 
Lodders et al.l ( l2009aU bl) also uses Si. For this reason we have plotted all of the data sets 
(Halley results, CI results and com et spectral rnqdelin g results) with both normalizations. 



Overall, the abundances derive d bvlLaw 



the solar photospheric values ( lAsplund 



er et al.l (Il989[ ) are close to or slightly smaller than 



20091 ). and the short and long data are also fairly 



consistent with each other for the heavier species. The data of I Jessberger et al.l (119881 ) and 
Jessberger fc Kisselll ( Il99ll ) exhibit anomalously high (above solar) values for Si and S when 
normalized to Mg. By comparison, normalizing to Si shows no excess over solar (except for 
Ca where the uncertainties are larger than most species). 



Which is to be preferred? The greater scatter in the data of IJessberger fc Kisselll (Il99ll ) 
may be due to the fact that normalizing to Mg versus Si will cause a scale change, as 
their Si/Mg ratio was ~2, and the uncertainties in converting the ion yields into atomic 
abundances were relatively uncertain. It may also be that the 74 sp ectra they us e d wer e 
not as representative of the Halley dust as the 433 spectra used by ILawler et al.l ( 119891 ). 
Certainly the confidence in the mean values of the abundances derived are likely to be more 



-36- 



representative of the true values for the larger data set, if they were otherwis e of co mparable 



quality. For these reasons we believe that the Halley data by iLawler et al.l ( 119891 ) and the 
normalization to Si are to be preferred. This choice doe not affect our main conclusions 
regarding the abundances of C and O, however. 
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Table 1. Spitzer IRS Observing Modules 



Module Abbreviation Wavelength Range (/im) Resolving Power {R = A/AA) 

Short Low 2"^^ Order SL2 5.2-8.7 60-127 

Short-Low P* Order SLl 7.4-14.5 60-120 

Long-Low 2"^^ Order LL2 14.0-21.3 57-126 

Long-Low P* Order LLl 19.5-38.0 58-112 

Short-High SH 9.9-19.6 ^ 600 

Long-High LH 18.7-37.2 ^ 600 



Table 2. Spitzer IRS Stare Exposure Times for Comet 73P/Schwassmann-Wachmann 3 
Slit Ramp Time (sec) Cycles Total Time (sec) Comments 



SL2 


240 


2 


480 


both nominal 


SLl 


60 


2 


120 


both nominal 


LL2 


120 


2 


240 


B saturated 16-19 /im 


LLl 


120 


2 


240 


both nominal 



-44- 



Table 3. Spitzer IRS Observations of Comet 73P/Schwassniann-Wachmann 3 



Object Mode UT Date UT Start Time Exposure (sec) Comments 



SW3-C 


SL2 


17 March 2006 


00:58 


SW3-C 


SLl 


17 March 2006 


01:17 


SW3-C 


LL2 


17 March 2006 


02:04 


SW3-C 


LLl 


17 March 2006 


02:11 


SW3-B 


SL2 


17 April 2006 


14:08 


SW3-B 


SLl 


17 April 2006 


14:29 


SW3-B 


LL2 


17 April 2006 


15:44 


SW3-B 


LLl 


17 April 2006 


15:24 



480 


nominal 


120 


nominal 


240 


nominal 


240 


nominal 


480 


nominal 


120 


nominal 


240 


saturated 16-19 fira 


240 


nominal 
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Table 4. NASA/IRTF SpeX Observations of Comet 73P/Schwassniann-Wachniann 3 - C 



Object 


Mode ^ 


UT Date 


UT Time 


Hour Angle 


Airmass 


SW3-C 


SXD 


20 April 2006 


10:56 - 11:47 


-1:01 - -0:21 


1.01 - 1.04 


SW3-C 


LXD 


20 April 2006 


09:47 - 10:42 


-2:11 - -1:02 


1.06 - 1.16 


HD 25152 


SXD 


22 August 2006 


14:25 - 14:34 


-1:55 - -1.46 


1.16 - 1.14 


HD 25152 


LXD 


22 August 2006 


14:11 - 14:20 


-2.09 - -2.00 


1.17 - 1.19 


^SXD = 


Short wavelength cross-dispersed echelle mo( 


le; LXD = Short wavelength 



cross-dispersed echelle mode 
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Table 5. NASA/IRTF BASS 3-13 /xm Observations of Comet 
73P/Schwassinann-Wachmann 3 



Object 



UT Date 



UT Time 



Airmass 



SW3-B 


18 May 2006 


15:33 - 


16:14 


1.13 - 1.06 


/SPeg 


18 May 2006 


16:20 - 


16:30 


1.06 - 1.05 


SW3-C 


18 May 2006 


16:35 - 


18:30 


1.05 - 1.03 


/3Peg 


18 May 2006 


18:37 - 


18:45 


1.04 


/3Peg 


19 May 2006 


15:30 - 


15:33 


1.15 


SW3-C 


19 May 2006 


15:43 - 


15:54 


1.15 - 1.18 


SW3-B 


19 May 2006 


16:13 - 


17:32 


1.01 - 1.08 


SW3-C 


19 May 2006 


17:39 - 


18:02 


1.02 


/SPeg 


19 May 2006 


18:07 - 


18:18 


1.02 


SW3-C 


19 May 2006 


18:24 - 


19:02 


1.04 - 1.07 
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Table 6. NASA/IRTF BASS CCD Hale-Bopp Blue Continuum (B^; 0.44 /im) & Red 
Continuum (R^; 0.71 /xm) Filter Observations of Comet 73P/Scliwassmann-Wachmann 3 



Object UT Date UT Time Airmass 



HD 164852 


19 May 2006 


14:30 - 


14:33 


1.13 


HD 170783 


19 May 2006 


14:40 - 


14:43 


1.13 


73P-C 


19 May 2006 


14:55 - 


15:10 


1.34 - 1.30 
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Table 7. Composition of the Best-Fit Model to the Spitzer IRS Spectrum of SW3-C 



Species 



Weighted " 
Surface Area 



Density 
(g cm"-^) 



M.W. 



^' moles 

(rcl.) 



Model Tmax" 
(°K) 



Model xl '' 
if not included 



Olivines 

Amorph Olivine (MgFeSi04) 
Forsteritc (Mg2Si04) 
Fayalitc (Fe2Si04) 

Pyroxenes 

Amorph Pyroxene (MgFeSi2 06) 
FerroSilite (Fe2Si2 06) 
Diopside (CaMgSi2 04) 
OrthoEnstatite (Mg2Si206) 

Phyllosilicates 

Smectite Nontronite 
Nao.33Fe2(Si,Al)40io{OH)2*3H20 

Carbonates 

Magnesitc (MgCOs) 
Siderite (FeCOs) 



0.00 


3.6 


172 


0.00 


225 


1.01 


0.13 


3.2 


140 


0.30 


225 


7.98 


0.02 


4.3 


204 


0.04 


225 


1.21 


0.34 


3.5 


232 


0.51 


225 


>100 


0.00 


4.0 


264 


0.00 


200 


1.01 


0.02 


3.3 


216 


0.03 


225 


1.28 


0.00 


3.2 


200 


0.00 


225 


1.01 



0.03 



2.3 



496 



0.01 



225 



1.74 



0.03 


3.1 


84 


0.11 


225 


1.41 


0.00 


3.9 


116 


0.00 


225 


1.01 



Metal Sulfides 

Niningeritc (MgsoFcyoS) 



0.06 



4.5 



78 0.35 



225 



3.43 



Water 

Water Ice {H2O) 
Water Gas (H2O) 

'Organics' 

Amorph Carbon (C) 
PAH {C10H14) 



0.15 


1.0 


18 


0.83 


200 


17.6 


0.02 


1.0 


18 


17 


200 


1.04 


0.53 


2.5 


12 


11 


280 


>100 


0.00 


1.0 


<178> 


<0.008 


NA 


1.01 



''Surface Area Weighting(i); A ± 25% 2-a error ellipse for the relative surface area of each species was determined by 
noting the fractional amount of change required to increase the reduced chi-squared value to just above the 95% confidence 
limit. 

^^moles ~Density(i)/Molecular Weight(i) 

"LTE at heliocentric distance of 1.47 AU = 233 K 

<*Total best fit model xl = 0-99; 95% confidence limit 



1.13 
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Table 8. Composition of the Best-Fit Model to the Spitzer IRS Spectrum of SW3-B 



Species 



Weighted " 
Surface Area 



Density 
(g cm"-^) 



M.W. 



^' moles 

(rcl.) 



Model Tmax" 
(°K) 



Model xl '' 
if not included 



Olivines 

Amorph Olivine (MgFeSi04) 
Forsteritc (Mg2Si04) 
Fayalitc (Fe2Si04) 

Pyroxenes 

Amorph Pyroxene (MgFeSi2 06) 
FerroSilite (Fe2Si2 06) 
Diopside (CaMgSi2 04) 
OrthoEnstatite (Mg2Si206) 

Phyllosilicates 

Smectite Nontronite 
Nao.33Fe2(Si,Al)40io{OH)2*3H20 

Carbonates 

Magnesitc (MgCOs) 
Siderite (FeCOs) 



0.04 


3.6 


172 


0.08 


245 


3.67 


0.09 


3.2 


140 


0.21 


245 


6.86 


0.00 


4.3 


204 


0.00 


245 


0.99 


0.28 


3.5 


232 


0.42 


245 


>100 


0.00 


4.0 


264 


0.00 


200 


0.99 


0.04 


3.3 


216 


0.06 


245 


2.50 


0.00 


3.2 


200 


0.00 


245 


0.99 



0.04 



3.2 



496 



0.00 



245 



0.99 



0.04 


3.1 


84 


0.15 


245 


1.85 


0.03 


3.9 


116 


0.10 


245 


1.65 



Metal Sulfides 

Niningeritc (MgsoFcyoS) 



0.01 



4.5 



78 



0.06 



245 



1.04 



Water 

Water Ice {H2O) 
Water Gas (H2O) 

'Organics' 

Amorph Carbon (C) 
PAH {C10H14) 



0.06 


1.0 


18 


0.33 


200 


3.50 


0.14 


1.0 


18 


120 


200 


1.93 


0.65 


2.5 


12 


14 


300 


>100 


0.04 


1.0 


<178> 


<0.023 


NA 


1.29 



''Surface Area Weighting(i); A ± 25% 2-a error ellipse for the relative surface area of each species was determined by 
noting the fractional amount of change required to increase the reduced chi-squared value to just above the 95% confidence 
limit. 

^^moles ~Density(i)/Molecular Weight(i) 

"LTE at heliocentric distance of 1.20 AU = 257 K 

<*Total best fit model xl = 0-99; 95% confidence limit 



1.13 
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SW3-B 



SW3-C 




SW3-B 



SW3-C 



Fig. 1. — Spitzer/IRS peak-up 16 /iin image maps of SW3-B (left) and SW3-C (right). 
Contours are plotted from 10, to 100 MJy sr~l (inclusive), in linear intervals of 18 MJy 
sr~^. North, east, the comet's projected heliocentric velocity vector (v), and the projected 
direction of the sun (©) are shown above each fragment. The image scale is indicated at the 
top, and is the same for both image maps. The morphology of the SW3-C coma appears 
similar to that of the majority of comets, while the SW3-B coma is highly elongated and 
diffuse, most likely due to the recent fragmentation event preceding the observations. See 
the electronic edition of the Journal for a color version of this figure. 
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Fig. 2. — Our computed slit-losses for our model fragment B surface brightness distribution, 
and those of a point source as provided by the Spitzer Science Center in the SPICE software's 
calibration files. Inset: Our model fragment B image convolved with a 16 /im peak-up array 
PSF. Contours mark the locations of the IRS 5-14 /im (SL) and 14-38 /im (LL) slits. 
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Fig. 3. — The same as Figure 2, except for fragment C The high-resolution 10-20 /xm (SH) 
module observation is also included. 
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Fig. 4. — The IRS spectrum of 73P/SW3-B. The lower (gray) hne is the spectrum normahzed 
to the flux extracted from the short wavelength slits. The upper (black) line is the IRS 
spectrum scaled by a fact or of 1.9 to niatch t he extracted flux from the 16 /im IRS imaging. 
The Gemini data i s from iHarker et al.l (J2010l ). The Spitzer IRAC and MIPS photometry is 
from lReach et al.l ( 120091 ). The photometry points have all been scaled to the same effective 
aperture size (p = 2900 km, or 5 peak-up array pixels), assuming a p~^ coma profile. 
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Fig. 5. — The IRS spectrum of 73P/SW3-C. The lower (gray) hne is the spectrum normahzed 
to the flux extracted from the short wavelength slits. The upper (black) line is the IRS 
spectrum scaled by a fact or of 2.4 to niatch t he extracted flux from the 16 /im IRS imaging. 
The Gemini data i s from iHarker et al.l (J2010l ). The Spitzer IRAC and MIPS photometry is 
from lReach et al.l ( 120091 ). The photometry points have all been scaled to the same effective 
aperture size {p = 5200 km, or 5 peak-up array pixels), assuming a p~^ coma profile. 
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Spitzer SW3-C Comet Spectrol Model 
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Fig. 6. — The spectral decomposition of SW3-C. In both this component and in B, the 
material of pyroxene stoichiometry is largely amorphous, while the olivine material is largely 
crystalline. The shortest IRS wavelengths are dominated by emission from amorphous car- 
bon. 
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Fig. 7. — The spectral decomposition of SW3-C, after removal of the emission by silicate 
materials. 
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Fig. 8. — The spectral decomposition of SW3-B. 
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Fig. 9. — The spectral decomposition of SW3-B, after removal of the emission by silicate 
materials. 
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Fig. 10. — The particle size distributions for the spectral models of SW3. 
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Fig. 11. — Mineral composition of components B and C The individual abundances for both 
comets are listed in Tables 7 & 8. 
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Fig. 12. — Water vapor bands in B and C compared. As discussed in the text, the strength 
of water vapor emission was highly time-dependent in both objects, so that any difference 
in strength was hkely simply a matter of the timing of the observations. 
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Fig. 13. — Atomic abundances of SW3 components B and C, compared to solar, normalized 
to unity for silicon. Also shown are the results for three other comets based on spectral mod- 
els, the CI carbonaceous chondrites, and IP /Halley (based on mass spectrometer results). 
The solar ph otospheric values are t aken from lAsplundl (120091 ). while the CI chondrite values 
are those of iLodders et al.l (l2009al Jbl). The uncertainties for the comet spectra models are 
two standard deviation values . 
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Fig. 14. — Atomic abundances of SW3 components B and C, compared to solar, normalized 
to unity for magnesium. 
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Fig. 15. — The relative timing of the mid-IR observations. The pre-perihe lion IRS and 
the B ASS data are discussed in this paper. The Michelle observations are in iHarker et al. 
(l2010i ). Post-perihelion data were also obtained with the IRS, but are of lower quality than 
the pre-perihelion data, and are not discussed in this paper. 
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Fig. 16. — SpeX observations of SW3-C. The data were obtained using the short wave- 
length cross-dispersed mode (0.8-2.4 Atm; SXD) and the long wavelength cross-dispersed 
mode (2.3-5.4 /im; LXD2.3) (JRayner et al.l 120031 ). Going from the shortest to the longest 
wavelengths, the spectrum transitions from being dominated by scattered solar radiation to 
thermal emission from the dust grains. 
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Fig. 17. — BASS observations of 73P-B. The observations have been broken down into 
individual sets that were interleaved with similar observations of SW3-C and flux calibration 
standard stars. The data on 19 May labeled "noscale" are an unweighted mean of all of 
the individual data sets during that time interval. Because errors in guiding, particularly 
during the daytime, tend to cause the throughput of the light to diminish (such errors never 
add light to the beam) the true signal may be somewhat larger than a simple average would 
indicate. Here "rescale" is an example of scaling the individual data sets with the highest 
one (integrating over the same wavelength range for each one). The difference between the 
two is a rough measure of the systematic uncertainty in the absolute flux density of the 
comet data, measured though the same fixed entrance aperture. 
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Fig. 18. — BASS observations of SW3-C. In addition to the changes seen between 18 and 
19 May, the data suggest significant intra-night changes. Note that twihght began at ~ 
1500 UT and sunrise at ~ 1600 UT, when guiding was done by centering the comet in the 
aperture to maximize the signal data sets. For component B all three data set form 18 May 
are identical, despite the second one being obtained in twilight and the third over 2 hours 
after sunrise. This suggests that the differences seen here in C are likely real, and not due 
to centering errors during daylight hours. 
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Fig. 19. — Components B and C, normalized to the same flux density at 12.5 /im. The 
relative strength of the silicate bands is consistent with the Spitzer IRS spectra. 
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Fig. 20. — Toy model of SW3-C using SpeX and BASS data to estimate the scattering 
contamination of the IRS data at 5 fira. The model was used to estimate the degree to 
which light scattered by the dust contaminate the shortest wavelengths of the IRS data used 
in the spectral modeling, and to help determine the grain albedos. 
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Fig. 21. — A{6)fp versus observation date for SW3-C. The 1'.'7 radius aperture data is shown 
as a +. The Re filter's A{6)f p extrapolated to a 4000 km aperture is shown as a x ( see fo r 
details). The squares are Johnson i?-band A{6)fp measurements from lBertini et al.l ( 120091 ). 
along with a qualitative trend line. 



